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An all-order resummation is performed for the effect of the running of the coupling a s 
in the zero recoil sum rule for the axial current and for the kinetic operator n 2 . The 
perturbative corrections to well-defined objects of OPE turn out to be very moderate. 
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The recent progress in the theoretical description of heavy flavors is based on 
application of the heavy quark expansion to QCD. Its main elements are nonrela- 
tivistic expansion and the Operator Product Expansion (OPE) allowing a treatment 
of the strong interaction domain in model-independent systematic way The precision 
achieved in certain cases is high and requires already a careful field-theoretic defini- 
tion of effective operators beyond a simple quantum-mechanical (QM) description. 
Such completely defined operators can be introduced in different ways, however, the 
common feature is their scale dependence. Likewise, the coefficient functions are al- 
ways /^-dependent. The general idea of separating two domains and applying different 
theoretical tools to them was formulated long ago by K. Wilson ]l| in the context 
of problems in statistical mechanics; in the modern language, applied to QCD it is 
similar to lattice gauge theories. The treatment of essentially Minkowskian quantities 
has, however, some peculiarities. They were first considered in || and then in 

In this paper I apply this technique to the zero recoil sum rules [§], to calculate 
one loop-corrections accounting completely for the effect of running of the strong 
coupling, which can be called an extended BLM || approximation, or, in short, merely 
BLM approximation. It is thus a direct BLM-generalization of the perturbative 
calculations of 0. In this way the BLM-improved perturbative evolution of the 
kinetic energy operator Q{iD) 2 Q = 7r 2 (/i) and the Wilson coefficient function 
for the zero recoil sum rule for axial current are obtained. 



1 Zero recoil sum rules 

The heavy quark sum rules and the method to calculate the relevant perturbative 
effects to them have been discussed in detail in 0; I will quote here only necessary 
equations. 

The first sum rule for the spatial component of the axial current at zero recoil 
through terms terms 1 / m 2 Q has the form 

/r M = (M AM ~/ M + A + it-) (i) 

6 m A c 4 \m c m b om c mb ) 

Here Iq A is the zeroth moment of the first structure function w± of a heavy hadron 
for the axial current: 

w A = 2hmh? t 
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= ^ rwf(e)de , e = M B -M D *-q Q . (2) 

Z7T JO 

Here and in what follows the ultraviolet cutoff in the moments is assumed to be 
introduced via 9{e — /i); of course, /i <C mq must hold. I use the standard notations 
for the expectation values 

^ = 2M^ {HQl ^ lHQh ^ = ^^(H Q \Q^G a PQ\H Q ) . (3) 

The sum rule and the lower bound for fi 2 was obtained using the sum rule for the 
vector structure function w\] for simplicity of notations below I consider the pseu- 
doscalar "weak" current J p = ci^^b and the corresponding structure function w p 
(the result for any similar current is the same): 

w ^j :<( ^, e ^ ) + ^^ym^M. (4) 

The Wilson coefficients £ depend on the concrete field-theoretic definition of fi 2 and 
Hq (in what follows terms ~ a sAq CD /mg are neglected and, therefore, n 2 G can be 
considered /i- independent, n 2 G ~ 3(M^» — M^)/4). At tree level one has 

u(ti = i , rV) = o . (5) 

The first - and so far the only - definition of the kinetic operator in field theory 
was given in 0: the expectation value was defined to have perturbatively = 
at any //. Thus, to calculate fi 2 in an external gluon field one must solve the Dirac 
(Pauli) equation of motion for nonrelativistic heavy quark, find the spectrum and 
integrate the spectral density induced by the current an from to fi. For an external 
soft field with frequencies much below /i this yields the classical value of (iD) 2 [0, 
therefore this is a proper definition of the operator n 2 in the quantum field theory. 
In general it clearly depends on /i. Given the definition, one is able to calculate its 
/x-dependence (mixing with unit operator) and the value of £4 without neglecting 
terms ~ fi 2 /rriq in the latter. 

2 Calculation 

The technical method was explained at length in & |3| ; it reduces to considering the 
OPE expansion in perturbation theory. Moreover, /Iq vanishes in perturbation theory 
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to leading order in 1/mg; it will be discarded below. In general, the kinetic operator 
is singled out by averaging operators over the hyperfine (Q-spin related) multiplet 
whereas the chromomagnetic operator is projected out. 

The new element here is the calculation of effects of running a s . In a one loop 
Euclidean calculation it is done by replacing a s by a s (k 2 ) where k is the gluon mo- 
mentum. For Minkowski quantities, in particular those that are not obtained by 
analytic continuation from Euclidean space, one should use the dispersive approach 
(it was used as long as 20 years ago and recently attracted renewed attention; for the 
extensive analysis and recent discussion see [§, |7j). It relies on using the dispersion 
relation for the dressed gluon propagator: 

r d\ 2 Ima s (-A 2 ) / k^k u \ 
J AT WT\ 2 \^- C ^) ■ (6) 

Then (5^ — ck^k y / k 2 ) / {k 2 + A 2 ) is a propagator of a gluon with mass-squared A 2 (the 
longitudinal component does not contribute in the one loop diagrams), and 

p = --lma s (-\ 2 ) (7) 

7T 

plays the role of the weight function for integrating over A 2 . More details are found 
in Q and papers mentioned there. 

As the first step a calculation of all quantities of interest in the one loop approx- 
imation with a non-zero gluon mass is needed. At A 2 = the calculations have been 
done in || and the modifications are straightforward. The perturbative spectral den- 
sities in the Ihs of the sum rules are given by the "elastic" peak of b — > c transition 

with e = and by the transition into the final state c+g(k) with e = u(k) = \J \ 2 + k 2 
starting e m i n = A . The effect of the quark recoil is of higher order in I/tuq and is 
discarded for l/m% terms. In the sum rule for p, 2 

1 rt 1 / 1 1 \ 2 

— / wf(e)de = p 2 M (8) 

2tt Jo 1 W \2m c 2mJ ^ VP; V ; 

the elastic transition is forbidden by parity. 

The inelastic transition amplitudes given by Figs. 1 are modified minimally com- 
pared to calculation of [0 but require certain care. One should distinguish between 
to and | A; | and use the proper two-body phase space |A;| doo; the scalar denominators 
of the propagators are ±2mQLj. Finally, the massive gluon has three polarization 



asjk 2 
k 2 



knkp 



k 2 
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propagating, instead of two for the massless one. The following trick simplifies the 
calculation. 

Since the tree gluon emission amplitude is transverse, the gluon propagator can 
arbitrarily be redefined in the following way: 



8 



(9) 



Choosing a = a = (k Q , — &)/2/cq = kj2k\ we arrive at the the tensor part of the 
propagator 



1 A 2 - - 

g»v = Sp, - ^(kftk v + k v k^ = e v± e l± + j2 n v n » > n=(Q,k/\k\) . (10) 

lKn i=l,2 







o 



In doing so we explicitly excluded the Coulomb quanta. The last term gives the 
contribution of the longitudinal polarization. In this form the rest of the calculation 
is simple. 

For the perturbative sum rule (H) we get 



— w p (lu; X)duj = 9{uj — A) 

Z7T 



2a s 



37r \m 



1 



1 



m b J 



(11) 



the effect of the gluon mass is only kinematic besides the longitudinal contribution 
A 2 /2a; 2 in the last bracket. 

As in the original analysis one can consider the sum rule for ^ k x ^ k currents: 



1 

27 



wX{e) de 



v 2 M - i*g(p) 
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1 1 

m c m b 



+ 2 



1 



3m 2 . 



• (12) 



The OPE states that all soft physics is given by the rhs ; the perturbative in- 
dependence of the moment is thus given by that of the operators appearing in it. 
Differentiating over ii we would get the perturbative spectral density, which thus 
must have the same functional form as eq. ([11]) with the proper dependence on quark 
masses. A direct computation yields 



2n 



w 1 (u; A) 



2a, 
9vr 



(lu-X) 



1 



m b 



-"-)' 

m c m b J 



1A 2, 
2 a; 2 



\k\ 



(13) 

the OPE, of course, is not violated in the one loop perturbative calculations.^] 

1 There were a few explicit checks || and even attempts to disprove |J OPE at this level in the 
recent years. 
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An (apparently) alternative definition of a4(//) was also given in p], [TXJ] using the 
third sum rule at v 7^ in the small velocity (SV) kinematics \v | -C 1: 

T^^W^ + ^W < 14 > 

3 Jo W ( e > v ) de 

where w is generated by an arbitrary heavy quark weak current having a non- 
vanishing nonrelativistic limit at v — 0. In this case the calculation is even simpler 
since the soft gluon radiation factorizes and for spacelike polarizations is merely given 
by the classical bremsstrahlung amplitude 

Ai oc g s ■ — x color factor . (15) 

UJ 

With e = uj in the SV limit one immediately arrives once again at the same one-loop 
contribution 

'f-i^t 1 ^)-* 1 - (16) 

in accord with the OPE. 

A similar calculation for the zero-recoil axial current sum rule (H) yields a more 
cumbersome result. The elastic peak is present and results in a contribution 

^<2(e) = 6(e) ■ (l + ^ro(A)) (17) 



where r (A) has been calculated in Jy], eq. (B.3). The calculation of the continuum 
part was detailed above and yields 



1 A (cont) / \ a s n , \ \ 

2^ crt W = ^-A) 



,1 1 2 \ A; 2 / 1 1 \ 2 A 2 fc 2 

2( — + — + 

b 





f 1 








- M 






m 6 y 



m 2 m? 3m c mb J uj 2 * ' ' ,2 ' ' 2 



1*1- 
(18) 



2.1 Resummation of perturbative corrections 

Integrals of the spectral densities eq. flTT|) and eqs. (ji7|) , QTH) are readily calculated: 



Am in 2 - A 2 ! 3 / 2 

2 \ pert/. .. \N _ /W..2 \2\\f J ' A J 



</4> pert (/i; A) = 6 (fx - A )— ^— . (19) 

5 



Inserting this result into the rhs of the sum rule (|l|) one gets 



£ A (W A) = 1+ 



a-* 



, r (A)+Ji — + — 



i2l _n + JjF 



-2A 2 log ■ 



"A 7 



A 



2(^_ 
3 



A 2 ) 3/2 



+ 



-> r (A| 2_ A 2 )5 /2 (// 2 _ A 2)3/2 



2W/i 2 - A 2 



#Gu 2 -A 2 ) . 



Using the explicit expression for r (A) |TTJ one finds 

3 (m b + m c m b 8\ . 
ro(A) = 7 log --A^log 



m Q \m c nib J 



(20) 
(21) 



£a(A) does not contain nonanalytic in A 2 terms through order l/mg which shows the 
absence of l/rrin infrared (IR) renormalon singularity in the BLM calculation. This, 
of course, is ensured by OPE since the infrared contribution below fi is peeled off 
from nonanalytic terms, on the other hand, can appear only in the infrared - 
elsewhere the propagator l/(/c 2 + A 2 ) is analytic in A 2 . The cancellation thus provides 
an independent check of the cumbersome expression for r (A). 

The absence of these nonanalytic terms ensures the absence of the corresponding 
IR renormalon in the resummation procedure; in other words, purely perturbatively 
the IR renormalons match in the sum rule. The concrete analysis therefore disagrees 
with the claims of of a discrepancy; basing on it the sum rules were declared 
erroneous in 0, |l!J. The criterion as applied in those papers does not make sense in 
general, and there is no point to address it here. 

One also notes the absence of nonanalytic terms in (/x 2 )i i 00 p, in spite of non-trivial 
operator mixing. It signifies the absence of the IR renormalon in Q(iD) 2 Q in the 
BLM approximation |5J (see also 0). The reason will be touched upon later. 

To calculate the BLM-type resummed result one must fix a particular form of the 
strong coupling according to one's preference. Then, for an observable A having the 
perturbative expansion through order a s 



A(X 2 



l + ^ l( A 2 

7T 



(22) 



2 The residual nonanalytic terms ~ A 3 /toq and smaller remain since the way to calculate £a 
is accurate up to such terms; they can be arbitrarily added or removed unless the operators of 
D = 6 and higher are incorporated. Their inclusion in the similar calculation would kill the next 
nonanalytic terms as well. 



6 



one has 

A™ = 1 + J A 1 (A 2 )p(A 2 ) . (23) 

For calculating £4 the exact infrared behavior of the coupling is not essential since 
the IR domain has been removed. 

The most popular choice for "all-order" resummation is the literal one- loop a s : 

^ = - ^ (24) 



l + £log£ Mog A 



V 



with 

P BLM (A 2 ) = \ [A V 5(A 2 +A y ) - 2 ; - - 1 . (25) 

Ay is Aqcd in the V'-scheme |@, A v = e 5/6 A^ D ~ 2.3A^ D ; here and throughout 

the paper a s is the V-scheme coupling: a s (k) = a^ IS (e _5/ ' 6 /c) neglecting higher-order 
effects in the ^-function. In this case one has 



A BLM = ^ + 4 




Ai(-A 2 v )- / ^A(A 2 



A 2 U 'log 2 ^ + 7^ 



(26) 



Using the explicit expression eq. ( p0|) the all-order BLM result for £a is obtained by 
simple integration. 

The /^-dependence of /x 2 and of the coefficient function £4 

d7r 2 (//) 

= 2c 7r (a a (/i)) fj, ■ I 



d/i 



2ae A (a s (//))/x (27) 



dyU 

in this approximation is given by integrating directly the continuum parts of w pert : 
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+2 (^-^) 2 ( 1+ ^-f)}FI" ( " 2 - A2) (28) 

Although the integrals here formally run over the "infrared domain" of A 2 as well, they 
actually do not include infrared effects due to specific constraints on the moments of 
p(A 2 ) |I3| with the weights in eqs. ( p^ ) analytic below fi 2 . 

In a similar way one could have tried to resum the perturbative series in /i 2 (/u) 
itself: 

A r rl\ 2 < n 2 — \ 2 W 2 

A^V) = (^l) pc M = ~ / ^ P BLM (A 2 ) (/i ^ 9tf - A 2 ) . (29) 

However, this procedure is physically senseless, no such object can be defined and 
any correction beyond BLM makes it ill-defined. It can be convenient, however, 
as an intermediate result to relate the particular field-theoretic definition of tt 2 to 
another one if the similar calculations in the latter are made. The absence of the IR 
renormalon in the one-loop approximation allows at least formal assigning a definite 
value to A^ 2 LM (/i). For this reason its anatomy will be addressed below. 



3 Analysis 

In Figs. 2 I draw the plots of A^ M (/i); Fig. 2a shows the dimensionless ratio r w 2 = 
Aj^| jM (/i)/A 2 / as a function of ///Ay; Fig. 2b illustrates the absolute value for three 
choices of \x. 
The value of 

VaM = (Urif 2 (30) 

is shown in Fig. 3; ?7a(/-0 must be used to calculate the physical formfactors in the 
OPE approach instead of T]a of HQET in the model estimates. Numerical evaluation 
was performed assuming m;, = 4.8 GeV and z = m c /mb = 0.3. 

For reasonable values of /1 and Ay natural for low-energy physics (the corre- 
sponding analysis |14j] suggests Ay ~ 300 MeV) the "perturbative" BLM-piece of 



/i 2 is clearly moderate, ~ 0.2 GeV 2 , and the latter does not depend too strongly on 
the renormalization point. Nevertheless, the BLM series, having a finite radius of 
convergence, is still divergent: the radius of convergence is given by 

f°M\ ^ = 1 ; (31) 

\ It J 7T 710 



3 The notations used here follow [la] where their meaning is clarified. 
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the series is sign-varying and asymptotically behaves like 



a s 7r 2 / im 3ir\ ( irb\ n ( a s (ii)\ n 

Even at Ay = 300 MeV the scale /x s h where the series starts to converge, a<j(/i s h) = 
4/6, constitutes 1.44 GeV. The first ten terms are written below: 

A BLM (/i) „ 4o^) . ^2 . N + 3 gg^. _ 2>51a; 2 _ 44>g;r 3 _ g ^4 + gg^S + g^G _ 

-25900x 7 - 40900a; 8 + 769400a; 9 + 162000a; 10 + ...1 ; x = ^1 (32) 

J 7T 

(this series is written in terms of the V^-scheme a s (//)). It is clear that the first BLM 
correction cannot represent numerically the complete result for any \x relevant for the 
heavy quark expansion in charm, as soon as the effects of the running of a s in this 
domain is taken into account. 

The analytic properties and the expansion are most conveniently obtained using 
the representation 

a blm/ \ 4a s (/i) 1 r dX 2 (y 2 - A 2 ) 3 / 2 a s (^ 2 ) 
A ^ W = —5— 7T~- f ~rr — r— r om...^ ( 33 ) 



3vr 2m IX 2 n i + ± as ^2 )hg 

with the contour shown in Fig. 4. 

The results for £4 are shown only in the limited range ji < 0.8 GeV; irrespectively 
of the actual onset of duality, values of \x above 0.7 GeV are not appropriate for OPE 
in charm: the power expansion runs in powers of /x/m c fe . With the running mass 
^c(a*) ^ 115 GeV one cannot sensibly adopt larger \i if the running of a s is included: 
in the BLM approximation the perturbative corrections are not small in this domain 
and the condition /i <C m c must be carefully respected. 

Assuming a s running one cannot allow larger values of Aqcd either; it is enough to 
recall that the "standard" value of Aq^ d ~ 250 MeV implies that the usual three-loop 
coupling hits the Landau singularity as early as at the gluon momentum ~ 900 MeV. 
These arguments having a more general practical relevance will be discussed in more 
detail in a separate publication. 

A clear illustration is merely the impact of 1/mq and higher IR renormalons in 

£4; as shown below already at Aqq D = 220 MeV one has 5 1 / m a (rf A ) ~ 0.05, thus 
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making the uncertainty irreducible without an account for 1/m 3 effect larger than 
the perturbative corrections themselves. 

The quantity C,a(^ 2 ', A) calculated through terms fj, 2 ' /m 2 still contains nonanalytic 
terms ~ A 3 , A 4 log A 2 etc.: 

ll + 5^ + 5^ 2 + ll^ 3 /A 2 \ 3/2 a. 

16 \ m l) n 

9 + 2z + 2z 2 + 2z 3 + 9z 4 A 4 A 2 a s , OA , 
Tc I lo S— •— • ( 34 ) 

The results shown in Fig. 3 correspond to a linear extrapolation of C,a{^ 2 ', A) in A 2 from 
A 2 = to —Ay which almost coincides with the principal value prescription; all arbi- 
trariness is an effect of Aq CD /m 3 and higher terms. The corresponding uncertainties 
(defined as the formal imaginary part) 



'1/m; 



(rf A ) ~ 0.0105 ( Av ) <Wfai) ^ 0.0015 ( Av ) 
KIAJ V300MeV/ Vrn Q v/AJ V300MeV/ 



(35) 

are still small at Ay = 300 MeV but too significant already at Ay = 500 MeV. 

Again, the radius of convergence of the BLM series for the Wilson coefficient 
(discarding in an arbitrary way l/mo and higher terms) is given by ay'(fx) = A/b ~ 
0.45 corresponding to a s (/i) — 0.29. 



4 Once more about the kinetic operator 

Certain confusion exists in the literature about the expectation value of the kinetic 
operator /i 2 . In the HQET a similar quantity is denoted — Ai; they coincide on the 
classical level but generally differ in the quantum field theory (QFT). However, a 
similar field-theoretic definition of — Ai has never been given. 

On one hand, — Ai is often defined as /i 2 (/i) from which some calculated per- 
turbative corrections are subtracted (different in different approaches) |i~6| ; on the 
other hand, the values deduced for such — Ai are also claimed for fi 2 [16|, [L7j - which 
certainly cannot be true simultaneously. 

The difference between possible definitions of the renormalized operator is seen 
already at order-a s perturbative corrections. They are given by two diagrams in 
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Fig. 5; both diverge quadratically and must be cut off. A direct computation shows [] 
that if one averages over directions of k a before integrating over k 2 , these diagrams 
exactly cancel each other, see eqs. (|A.1| ) - (|A.2| ). Therefore, if a cutoff is introduced 



by any function depending on k 2 , the one- loop renormalization vanishes. Clearly, one 
does not get renormalization in BLM resummation either. As a matter of fact, in the 
Abelian theory this holds up to the order a 3 , and without light flavors (m q ^ fi) in 
all orders in a (see Appendix). In the non- Abelian theory, on the other hand, such 
an "isotropic" averaging merely cannot be defined in a gauge-invariant way beyond 
one loop. 

The operator 7f 2 , however, is defined differently in |fjj. It is easy to see that at 
the one loop level it corresponds to introducing 6(fi — \k |) with no constraints on ko- 
In such a case the renormalization is obviously present, and is given by c n = 4a s /37r 

§■ 

The 'Lorentz-invariant' cutoff at k 2 = fi 2 at first sight seems more convenient. 
However, the advantage is purely technical and applies exclusively to the one-loop 
perturbative calculations.^ It is possible to work with such definitions in Euclidean 
theories and statistical mechanics. 

On the other hand, theory of heavy quarks is essentially Minkowskian, and a 
'Lorentz-invariant' cutoff generally does not allow to formulate theory on the nonper- 
turbative level. Imposing constraints on the timelike components of momenta leads 
to non-locality in time, and the standard Hamiltonian approach is not applicable. 
For example, a heavy hadron state \Hq) over which one wants to average Q(iD ) 2 Q 
is not completely defined. Using another language, the Euclidean amplitudes either 
cannot be continued to the physical domain, or possess wrong analytic properties 
there. These deep theoretical questions will be discussed in subsequent publications. 

Practically, as soon as a proper particular definition of an operator is given, one 



4 This simple calculation has been checked in discussions with M. Shifman and A. Vainshtein in 
late 1994 to clarify the observation || about the disappearance of the renormalon in the kinetic 
operator in one loop. 

5 It is worth noting that both ways to introduce a cutoff are equally Lorentz-covariant: #(/i 2 — 
k 2 ) = 9(n 2 + k 2 — (vk) 2 ) where v is velocity of the heavy hadron. The difference is rather the 
explicit dependence on the hadron state (its velocity). However, this distinction does not have 
a general meaning since the kinetic operator itself Q {—D 2 + (vD) 2 } Q is, in this respect, not 
Lorentz-invariant either depending explicitly on the velocity v a . Alternatively, the nonrclativistic 
heavy quark field Q is not an invariant since is constructed from the QCD field Q in the non-invariant 
way Q = e lmQ ( vx ^Q and, thus, evolves according to a non- invariant equation of motion i vD Q = 
involving explicit vector v a . In the non- Abelian theory the attempt to define a u-independent cutoff 
on the gluon momenta is necessary a gauge-dependent procedure. 
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can relate its expectation value to that in other schemes. For example, in routine 
calculations in the QCD sum rules the cutoff is extrapolated to zero in the one- 
loop approximation. Such a value of — Ai is less than fi 2 by some small amount 
calculated in |2| , therefore, strictly speaking, in a number of practical applications /i 2 
must be taken even larger than estimates quoted in JT3]. On the other hand, for purely 
technical reason a definition with the covariant cutoff was applied in the calculation 



of the semileptonic widths in 0, therefore the literal numbers of [[18] were used 
there, which are somewhat smaller than jj^ discussed in || and in the present paper. 
The numerical analysis above easily allows one to use a more convenient calculation 
scheme in each particular case. It is worth emphasizing again || that subtracting 
the perturbative piece from the renormalized operators does not make sense by itself 
when nonperturbative corrections are addressed, and the result is usable only as long 
as the amount subtracted from the expectation value of properly defined operator is 
explicitly specified. 

The advantage of the "physical" definition of fi 2 is that it allows to derive an 
important lower bound 

A > /4 (36) 

and still leads to moderate - and well-controlled - perturbative corrections in observ- 
able quantities. In any case, so far no sensible alternative definition has been given 
in the literature; in principle, of course, it can be done. 

It is important to note that the relation between properly given QFT definitions 
of an operator is done perturbatively and excludes the IR domain. For example, if 
there is no renormalon in the perturbative series for it in a certain approximation in 
one scheme, the same holds in any other scheme as well. This explains the absence of 
the renormalon in the BLM series for /^(/i) noted in Sect. 2.1: with the cutoff over 
k 2 all terms merely vanish. 

It does not contradict the presence of mixing: 

d 2 = <Ma.s(/i)) • (37) 

Even in the large-n^ approximation whole series in a s (fj,), see eq. (EBf) . It is 

easy to understand this fact: calculating (j> 2 ([j>) perturbatively one integrates over the 
'cylinder' domain 

— oo < k < oo , \k\ < fj, . 
On the other hand, the integral over the spherical domain k 2 < fj 2 vanishes. For this 
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reason the integral can be reduced to the domain 



ko > \l [i 2 — k 2 } \k\< n 



(see Fig. 6) where k 2 > fi 2 always holds and the coupling remains small. 

It means that in eq. fl37|) the coefficient c n is given by an effective coupling that 
never grows to infinity at finite /i in the BLM approximation. These subtleties were 
missed in [19| where mixing eq. (|37|) was equated to the presence of renormalons. 

The situation is completely different in the case of the pole mass: the leading IR 
contribution bvrtq ~ [i comes from k ^ — ^ ~ [20|, therefore k 2 = —k 2 and 



the corresponding coupling is always a s {jj) in the BLM approximation. 



5 Conclusions 

All-order BLM effects are calculated for the zero recoil sum rules within proper OPE 
approach. Contrary to existing claims they have quite moderate impact. In particu- 
lar, for the available field-theoretic definition of the kinetic operator it 2 its expectation 
value obeys the inequality [|], 

P 2 M > /4 

for any normalization point fi; the normalization-point dependence is very moderate. 



The conclusion differs from a recent paper |16[ where the first nontrivial BLM 
terms 0(ba 2 ) were addressed; those represent the first term of the expansion of 
the full function calculated in the present paper. It is seen that the series for this 
functions are still divergent and the calculated terms already grow in magnitude; in 
this case the estimates based on the first term 0(ba 2 ) are numerically misleading. 
The analysis shows that for the scale governed by charm effects of running of the 
coupling - if addressed at all f] - must include the whole resummation, made in the 
OPE-consistent way. 

Moreover, regarding the perturbative corrections to the sum rule for the zero- 
recoil formfactor F D * (the axial-current sum rule ([[])), the calculation of 0{ba 2 {ji)) 
terms in fl6j missed the similar contribution from r) 2 A \ considering that calculation 
as the (second-order) perturbative correction to the sum rules is thus erroneous. It 

6 The phenomenological relevance of such theoretical improvement is far from obvious; it will be 
discussed elsewhere. 
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is only the whole set of corrections that eliminates the infrared domain from the 
perturbative factor. As a result of this omission the major part of the corrections 
came from the domain near and below the Landau singularity and therefore seems 
irrelevant. 

It is worth mentioning that the smearing procedure considered in [[UJ practically 
does not improve the convergence: the critical value of a s for the weight fi 2n / (fi 2n + 
e 2n ) decreases by only a factor of 1 — l/(2n) for the price of increasing the actual 
cutoff scale while the room for it is severely limited by the small value of m c . 

The value of rjA(/j) = (£a(/-0) 1//2 & t a representative choice // ~ 0.5 GeV, Aqq D = 
300 MeV is 

T]a{ij) = 1 tree level 

t)a{ij) — 0.975 one loop 
T] A (n) ~ 0.99 all-order BLM 

Clearly, the perturbative corrections including BLM effects are very moderate and 
differ minimally from the estimate 0.98 used in the original analysis [|], |2|. The 
difference falls well below the effect of 1/mj? corrections not addressed so far. 

The nonperturbative corrections, on the other hand, do have a grave numerical 
impact on the model analysis [JT^, |22| (see also |J) where it was postulated to use 



ill-defined fj A H( ^ ET \= ?7a(0) ) according to the routine practice of HQET P3"] . In all 
later publications the perturbative factor was used about 0.95 which is not supported 
by the analysis. Moreover, such an approach results in double-counting of the soft 
domain contributions. 
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Appendix 

At one loop the renormalization of /x 2 is given by Figs. 5. The dark box is the operator 



14 



Q(iD) 2 Q (0) . The corresponding Feynman integrals are 

f d 4 k 1 1 r2 T r d 4 k 3 



(A.l) 




(A.2) 



Thus, with such a cutoff I a + h = 0. If, to preserve the QM interpretation, one first 
integrates over ko from — oo to oo , one obviously has 3I a = —If, . 

Let us consider higher orders in QED without light flavors. The Lagrangian is 



where the (Euclidean) Lagrangian of the Abelian gauge theory with the isotropic 
cutoff at k 2 = /i 2 (/i <C itlq) is defined as 



The dynamic degrees of freedom are A^{k) with k 2 < ji 2 . Note that, formally, the 
gauge invariance is lost due to the constraint on A^; the current is still conserved, 
however. In this theory ^2,3 are sterile; only Aq interacts with Q. Thus in higher 
orders one has only the Coulomb-exchange dressing of two one-loop diagrams Figs. 5, 
illustrated in Figs. 7a and 7b, respectively. It will be shown below that both are 
renormalized multiplicatively by the same factor (which, as a matter of fact, equals 
unity) and thus the overall renormalization of the kinetic operator vanishes if the 
renormalization scheme is such that two first-order graphs in Fig. 5 cancel. 

This follows from the exact factorization in the emission of the Coulomb quanta. 
Since the gauge field self-interaction is absent (generally, it is a cutoff-dependent 
fact), the perturbation theory is fully time-ordered in the heavy quark limit. It is 
convenient to classify the perturbative diagrams by the number n of Coulomb quanta 
spotted at t — 0. Consider first the tree-level emission of n quanta from the initial 
state (terms e k for the amplitude with k = n). The corresponding sum over possible 
permutations is given by 




(A.3) 



d A k 



A„(k)Al(k) (k 2 5^ - kpK) 9(^i 2 - k 2 ) + gauge fixing terms. (A.4) 



E 



1 



1 



1 



n 



1 



(A.5) 



permutations 



OOi UJi +U 2 U1+UJ2 + UJ 3 
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All virtual emissions cancel between the wavefunction and vertex renormalizations 
due to the QED Ward identity [Z\ = Z 2 ), and possible formfactors are absent in the 
limit [ijviiQ — > 0. Thus, in this static limit the emission amplitude for n photons is 
exactly e n T7 ~r. to all orders. This property is nothing but the quantum decomposition 
of the classical Coulomb field of a static source. 

The Feynman integral for the sum of the diagrams Fig. 7a with n Coulomb quanta 
at t = is then given by 

J(n) = (47ra)" i g d% 



n\ J f- = \ (2n)H kfu, 




(4na) n r d 4 fc w k 2 r y} _J__ 
" (n-1)! J (2K)Hk*u* J i\ (2n)Hkfuf 1 ' ) 

whereas, with k n being the momentum of the spacelike-polarized photon in the blob, 
the expression for the diagrams Fig. 7b with n — 1 Coulomb quanta is clearly 

r(B) _ (47m)" f d*k^3^ f g _^L^_ . (A7) 



(n-1)! J (2rr)Hkl J £J 

Thus the diagrams Fig. 7a and 7b cancel order by order. Note that the common fac- 
tors in eqs. ( |A.6| ) and ( A.7|) are merely the vertex correction Z\ for the static charged 



particle to the given order, which coincides with the wavefunction renormalization 
Z 2 'i it cancels the diagrams with n = 0. Therefore the one- loop value of ^ simply is 
not renormalized by higher orders if the regularization yields a non-vanishing value 
in order a. This is also rather obvious starting from the definition eq. fll4]) . 

Since the result above does not depend on a particular form of the gauge propa- 
gator and only requires it to be an invariant function of k 2 , the insertion of fermion 
loops does not destroy the cancellation as long as effective multi-photon vertices do 
not appear. Due to Furry's theorem it occurs only at order a 3 when light fermions 
are present; examples are shown in Figs. 7c, 7d. 

A similar cancellation is not expected to hold in the non-Abelian case already 
at two loops [Q] where the difference with the Abelian case first emerges. However, 
the isotropic cutoff by itself cannot be defined in a gauge-invariant way; for example, 
the simplest Pauli-Villars regularization violates color current conservation, and the 
renormalization, instead, would depend on the gauge. There is a more complicated 
general method to introduce a similar regularization in the invariant way |24j]; the 
presence of mixing there is, however, the most general feature in two loops. 
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Figure Captions 



Fig. 1: Perturbative diagrams leading to nontrivial inelastic structure functions. 

Fig. 2: The perturbative part of the kinetic operator in the resummed BLM approx- 
imation, A^ 2 LM (/i). 

a) : The dimensionless ratio r w 2 = A^ M (yu)/ A v as a function of the dimensionless 
scale parameter ///Ay. 

b) : The absolute value of A^ 2 LM (/t) as a function of Aq2 d f° r A 4 — 0.5 GeV, \i = 
0.75 GeV and n = 1 GeV. 

Fig. 3: The value of the Wilson coefficient r} A {n) = {U{v)) 1/2 for A§2 D = 200 MeV, 
Aqq D = 300 MeV and Ajjgo = 400 MeV. The value of r] A {fi) must be used in the 
QCD-based calculations of the exclusive zero-recoil B — > D* formfactor when non- 
perturbative effects are addressed. The shaded bars show the purely perturbative 
uncertainty irreducible without the OPE account for \jw? c and l/m* effects; they 
thus represent the lower bound for the corresponding actual nonperturbative correc- 
tions. 

Fig. 4: The contour of integration over A 2 allowing a straightforward perturbative 
expansion of A^ 2 LM (/i) • 

Fig. 5: Feynman diagrams contributing to the one- loop renormalization of the op- 
erator Q (iD ) 2 Q . Dashed line is the gauge boson (gluon); dark box represents the 
operator. In the diagram a) only the Coulomb quanta (A ) propagation contributes 
whereas in b) the gluon is spacelike. 

Fig. 6: The domain of integration over the gluon momentum k; the cylinder —\i < 
\k\ < n gives the renormalization of /x 2 in one loop. The integral over the shaded 
disk (sphere) vanishes. 

Fig. 7: Higher-order renormalization of the kinetic operator in the Abelian theory, 
a, b: without light flavors (n — 2). 

c, d: presence of light flavors induces photon self-interaction in the order a 2 . 
Dashed photon lines denote Coulomb quanta, dotted lines are for the spacelike po- 
larizations. Dark box is Q (iD ) 2 Q(0) . Shaded blobs show the possibility to twist 
photon lines. 
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